Crossover from 2D to 3D in a weakly interacting Fermi gas 
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We have studied the transition from two to three dimensions in a low temperature weakly inter- 
acting ®Li Fermi gas. Below a critical atom number, N2D, only the lowest transverse vibrational 
state of a highly anisotropic oblate trapping potential is occupied and the gas is two-dimensional. 
Above N2D the Fermi gas enters the quasi-2D regime where shell structure associated with the filling 
of individual transverse oscillator states is apparent. This dimensional crossover is demonstrated 
through measurements of the cloud size and aspect ratio versus atom number. 
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Dimensional crossovers provide a bridge between seem- 
ingly disparate behaviors in physical systems. For exam- 
ple, Bose-Einstein condensation (BEC) and superfluidity 
generally coexist in three dimensions (3D) but, in 2D, 
Berezinskii-Kosterlitz-Thouless (BKT) superfluidity can 
occur without BEC [TJ|5]. Ultracold atomic gases offer 
unprecedented opportunities to explore such crossovers 
where the dimensionality and interactions can be pre- 
cisely controlled [3]. The evolution from a BEC to a 
BCS (Bardeen-Cooper-Schrieffer) superfluid at a Fesh- 
bach resonance is now a central topic in 3D Fermi gas re- 
search PHTU] . Restricting these gases to 2D dramatically 
alters pairing and superfluidity |11H19| and may offer in- 
sight into itinerant ferromagnetism |20l ET] and Fulde- 
Ferrell-Larkin-Ovchinikov superfluidity |22ff24] , also seen 
in quasi-2D heavy fermion superconductors [25] . In the 
crossover region, the evolution from a BCS to BKT su- 
perfluid can also be investigated. 

Dimensional crossovers are characterized by access to 
new degrees of freedom. As new quantum states become 
accessible, their discrete energies can be immediately ap- 
parent. In quasi-2D Fermi gases, where the transverse 
confinement energy is on the scale of the Fermi energy, 
occupation of new transverse states gives rise to shell 
structure, leading to steps in the density profile, chemi- 
cal potential and specific heat with increasing system size 
as a consequence of Pauli exclusion [26 - 28J . 

In this letter we study the crossover from 2D to 3D in 
a two-component optically trapped ^Li Fermi gas. The 
scaling of the cloud width with atom number in both the 
tight and weakly confined directions changes dramati- 
cally through the crossover. We also see clear indications 
of shell structure associated with the filling of discrete 
energy levels, a feature which is most pronounced in the 
aspect ratio of the cloud. The experimentally observed 
structure agrees well with theoretical predictions for a 
weakly interacting Fermi gas. 

Recent experiments have investigated 2D phenomena 
in Fermi gases prepared in one-dimensional optical lat- 



tices |29ff55] . Lattice-based experiments produce mul- 
tiple 2D clouds that can be imaged simultaneously of- 
fering good measurement signal-to-noise. However, the 
cloud size in the tightly confined direction is typically 
sub-micron, well below the resolution of nearly all imag- 
ing systems, so the properties of individual clouds cannot 
be easily measured. We overcome this limitation by cre- 
ating single 2D clouds, that can be expanded for imaging. 

Achieving the 2D regime in a Fermi gas requires that 
the Fermi energy, Ep, and temperature, T, are suffi- 
ciently low that excitations in one dimension {z) are for- 
bidden. This can be understood by considering an ideal 
Fermi gas conflned in an harmonic trapping potential 

V{x, y, z) = ]^m{ujlx^ + uj^y"^ + ^lz^)i (1) 

where m is the mass of the atoms and uJr,z are the trap- 
ping frequencies in the radial (r) and transverse (z) di- 
rections, respectively. For simplicity, we consider the ra- 
dially symmetric case — cOy — w^) and an oblate 
geometry where uJz 3> uJr- lu 2D the Fermi energy is 
given by Ef,2D — V 2Nhujr where N is the number of 
atoms in each spin state. In 3D the Fermi energy is 
Ef,3D ~ {6Ny/^hu} where w = (wj.Wj^w^)^/'^. 

At zero temperature, the Fermi radii are set by 
the Fermi energy. In 3D this gives Rp — 
(48iV) hui I (moj'^. ) , where = x,y,z. In contrast, 
when the gas is two-dimensional (Ep <^ hujz), the width 
in the tightly confined direction is set by the size of the 
harmonic oscillator ground state = h / [muj z) ■ In the 
radial directions, however, the Fermi radius is 

Rp^r = mfUf^. (2) 

V mujr 

This A^^/^ growth is more rapid than in 3D as two phase- 
space degrees of freedom {z^pz) are no longer accessible. 
As Ep depends on the atom number we define a critical 
number N2D below which atoms populate only the low- 
est transverse vibrational state at T = and the gas is 
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2D. This is equal to the number of single particle states 
with energy less than the lowest state with one trans- 
verse excitation. Labeling harmonic oscillator states by 
the vibrational quantum numbers rij. and and defin- 
ing the trap aspect ratio as A = uj^l'^r we can count 
the number of (n^ — 0) states with energy less than the 
{rir = 0, riz = 1) state. These states satisfy rir < X so, 
including degeneracies, N2D is given by 
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In the experiments that follow, we work in a trap with 
A « 60 corresponding to N2D ~ 1800. 

When only a few transverse states are occupied the 
gas is quasi-2D and the number Nn at which the n^^ 
transverse state begins to fill is given by 



Nn = yn(n+ 1) 
4 
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At these points we expect steps in the scaling of the cloud 
size with N , giving rise to shell structure. 

We quantify the cloud size using root mean square 
(rms) radii — \/ {ff) which provide a model indepen- 
dent width measure applicable at both zero and finite 
temperatures. To calculate theoretically, we first find 
the 3D density profile n{x,y,z), integrate this over two 
dimensions to obtain a line profile n(ri) and then evalu- 
ate the second moment (rf) = / n{ri)rfdri/ J n{ri)dri = 
a^.. For an ideal Fermi gas, n{x,y,z) is found by sum- 
ming the squared wavefunctions of the individual oscil- 
lator states. However, as our experiments are performed 
at a magnetic field of 992 G where the s-wave scattering 
length is — — 43G0ao {qq is the Bohr radius), we have 
also calculated n{x,y,z) by numerically solving the 3D 
Hartree-Fock mean-field equations for our oblate trap- 
ping potential. In the weakly confined radial directions, 
we use a local density approximation assuming a slowly 
varying density profile as a function of r. At each r, we 
then solve the Schrodinger equation for the transverse 
direction to obtain the full density profile n{x, y, z). 

To create a 2D Fermi gas, we begin with a cloud of 
approximately N = 10^ ^Li atoms in each of the lowest 
two spin states \F — 1/2, mp = ±1/2) in a far detuned 
optical dipole trap. The cloud is evaporatively cooled to 
a temperature T 0.1 Tp at a magnetic field 834 G ^4J, 
at the centre of the Feshbach resonance, where elastic 
collisions are unitarity- limited. To vary the final atom 
number we continue the evaporation by further lowering 
the dipole trap power so that atoms are spilled while 
maintaining the cloud at the lowest possible temperature. 

Next we adiabatically ramp on the 2D optical trap, 
formed by a cylindrically focussed Gaussian beam prop- 
agating along the y-direction with 1/e^ waists of Wz = 
8 /im and Wx = 400 fim in 200 ms. Similar configurations 
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Figure 1. One-dimensional density profiles for clouds at 992 G 
as a function of atom number in the tight (a) and weakly 
confined (b) directions after 500 /is expansion. Solid black 
lines are contours of equal density and the dashed red lines 
show the predicted rms cloud width for an ideal Fermi gas. 
The color bar indicates line density in units of atoms//im. 



have been used to study 2D Bose gases |351 [55] . Once 
this trap is fully on, the first beam used for evaporation is 
ramped down in 200 ms leaving the atoms in the 2D trap. 
Confinement in the y-direction is achieved by the short 
Rayleigh length associated with the 8 /im waist. Addi- 
tional confinement in the x and y-directions is provided 
by the curvature of the magnetic field. The 2D trapping 
frequencies are utz /2Tr = 2800 Hz and uj^I^-k = 47 Hz 
(cjx ^ iOy = ujr) giving an aspect ratio of approximately 
60. Finally we adiabatically ramp the magnetic field to 
992 G where the cloud is imaged. Precise temperatures 
in the 2D trap are difficult to ascertain due to the lack of 
an analytic model for interacting quasi- 2D gases. How- 
ever, we expect the temperature to be below 0.1 Tp due 
to the adiabatic field sweep and deep evaporation used 
to prepare the low atom number clouds. 

To demonstrate the crossover from 2D to 3D we mea- 
sure the cloud radii in the tight and weakly confined di- 
rections as a function of iV. In the ^-direction the cloud 
width in trap (a^ 770 nm) is much smaller than the 
resolution of our imaging system, so a short time of flight 
(500 /ts) is used before imaging. This time is long com- 
pared to the inverse trapping frequency in the z-direction 
{l/cjz — 57 /is) but short compared to l/wj, (3.4 ms). The 
cloud distribution in the radial direction will therefore be 
equivalent to the in-trap distribution. 

The imaging beam propagates roughly along the (ra- 
dial) a;-direction, so we obtain two-dimensional density 
profiles n{z,y). Absorption images are processed using 
a fringe removal algorithm |37] to optimize image qual- 
ity. From these images we generate line profiles n{z) and 
nly) by integrating over the other dimension. Figure 1 
shows the measured line-densities as a function of atom 
number for gases at 992 G (og = — 4300ao) in both the 
tight and weakly confined directions. The plots consist 
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of approximately 50 profiles; each profile is the average 
of several images binned according to N. The dashed 
red lines show scaled predictions for the zero tempera- 
ture rms radius for an ideal Fermi gas and the solid lines 
are contours of equal density. 

In order to compare these data with theory, we evalu- 
ate the rms widths from the profiles in Fig. 1. Figure 2 
(main panel) shows the rms cloud width in the j/-direction 
(cTy) versus atom number. Data points are the experi- 
mental measurements and the solid and dashed lines are 
the predictions for a weakly interacting and ideal Fermi 
gas, respectively. Each theory curve has been scaled by 
1.19 (1.11) for the weakly interacting (ideal) gas to lie 
on top of the experimental data at low TV, where interac- 
tions are least significant the models are most accurate. 
These scalings account for factors such as finite imaging 
resolution, nonzero cloud temperature and a slight ellip- 
ticity in the radial confinement (ux/i^y ~ 1-1), but do 
not affect the power-law dependence of the width. The 
data closely follow the predicted growth rate for a weakly 
interacting gas over the full range of atom numbers with 
an elbow around N = 2000 (~ N2D ) corresponding to 
the transition between 2D and quasi-2D regimes. The 
data deviate below the ideal gas prediction at high N 
where interactions become more important. As the atom 
number and density increases, the attractive interactions 
become more significant and slow the growth rate of the 
cloud. For the lowest numbers (iV « 800) the interac- 
tion parameter 1/ikpas) = —2.3, where kp is the Fermi 
wave- vector. For N k, IQ^ the gas is approaching the 
strongly interacting regime, l/(fcf^a<;) = —1.0, and our 
weakly interacting theory will begin to break down. Be- 
low N2D — 1800, a least squares fit to the experimental 
data shows that the width grows as 7V"-^®='=°-°^, in agree- 
ment with the N^-"^^ scaling predicted by Eq. (2). Above 
N2D1 the fitted dependence is TV'^ -'^^^^'^ "^'', in reasonable 
agreement with the interacting gas prediction of 7V''.i46 
but well below the ideal gas prediction of N^-^'^'^. 

The inset of Fig. 2 shows the cloud width in the z- 
direction. The theoretical widths here have been scaled 
by b^Jl -\- uj^t^ to account for the t = 500 fis expansion. 
The value of b which provides the best agreement with 
the low N data is 1.24 (1.20) for the weakly interacting 
(ideal gas) theory. These values are slightly larger than 
those used in the radial direction suggesting the expan- 
sion may not be simply ballistic. Below N2D the trans- 
verse width is roughly constant scaling as 7V*' °^*° '^"' in 
reasonable agreement with the predictions of iV^*^ "^ for 
an interacting gas and for an ideal gas. Finite tem- 
perature may lead to occupation of the (n^ = 1) state 
for N < N2D which would increase the measured expo- 
nent. For N > N2D the width scales as 7VO-i24±o.oo4^ 
compared to iV''-^°° and TV"-^^"^ for the interacting and 
ideal gases, respectively. These exponents measured in 
the z-direction may be infiuenced by the cloud expansion 
which could have a non-trivial dependence on N |38l [55] 
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Figure 2. Measured rms cloud widths (data points) and the- 
oretical predictions for a weakly interacting (solid lines) and 
ideal (dashed lines) Fermi gas. Radial cloud width is plotted 
in the main panel and the inset shows the transverse cloud 
width after 500 /is expansion. Theoretical curves are scaled 
as described in the text. 



as 1/ {kpUs) varies by more than a factor of 2 for the range 
of atom numbers considered here. We note that the cal- 
culated exponent for the non-interacting gas is below the 
true 3D value of 1/6 as the data lie in the quasi-2D regime 
where shell structure infiuences the widths. 

To further investigate the shell structure we now focus 
on the cloud aspect ratio, k — Uzjoy. Shell structure 
will be more prominent in measurements of n, than in the 
individual widths as the filling of new transverse shells 
begins in states with low radial quantum numbers; hence 
an increase in the transverse cloud size will correlate with 
a decrease in the growth rate of the radial size. In the 
3D limit, the cloud aspect ratio would be constant but 
in 2D and quasi- 2D k will show a strong dependence on 
the atom number. Additionally, as n, is given by the 
ratio of two measured quantities, certain experimental 
systematics (e.g. finite imaging resolution and shot to 
shot temperature variations) will be reduced. 

In Fig. 3 (main panel) we plot the aspect ratio of the 
cloud along with theoretical predictions for the weakly in- 
teracting (solid line) and ideal (dashed line) Fermi gases. 
The agreement between theory and experiment is very 
good and the large change in aspect ratio with iV clearly 
demonstrates departure from 3D behavior. Only at high 
atom numbers does the aspect ratio level off indicating 
broad coverage of the 2D and quasi-2D regimes. The ar- 
rows indicate the calculated atom numbers at which new 
transverse states become accessible, Eq. (4). The inter- 
acting and ideal gas predictions are very similar, empha- 
sizing the robustness of aspect ratio measurements for 
identifying the dimensional crossover. 

At low A'', the aspect ratio decays steadily with in- 
creasing atom number, before a step occurs at TSS2D cor- 
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Figure 3. Cloud aspect ratio in the 2D-3D crossover (points) 
and theoretical predictions for a weakly interacting (solid line) 
and ideal (dashed line) Fermi gas. Arrows indicate predicted 
atom numbers at which new transverse states become acces- 
sible. Inset: Gradient of the experimental and theoretical 
aspect ratios showing the signature of shell-filling. 



responding to occupation of the first transverse excited 
state. The inset of Fig. 3 shows the gradient dn/dN 
of the theoretical and experimental aspect ratios which 
show the signatures of shell structure more clearly. The 
experimental gradients were evaluated numerically and 
smoothed with a five-point moving average. Both the 
measured aspect ratio and gradient closely follow the the- 
oretical predictions, with indications of shell structure 
present in the experimental data for N < 10, 000 corre- 
sponding to occupation of the ground and first transverse 
excited states. The location of the first two arrows in the 
plots agrees well with the position of the steps. For larger 
N the data lie close to the theoretical line but the shell 
structure is unresolved. 

This work provides the first quantitative study of the 
transition from 2D to quasi-2D and 3D in a weakly in- 
teracting Fermi gas. At low atom numbers, shell struc- 
ture, associated with the filling of individual transverse 
oscillator states, becomes apparent. Our data were ob- 
tained away from any confinement induced scattering res- 
onances {qs < ttz); however, these could easily be ac- 
cessed closer to the Feshbach resonance [l2l [3l] 140] . This 
work opens the way to investigations of the phase dia- 
gram of 2D and quasi- 2D Fermi gases as a function of 
the scattering length and temperature and could help 
elucidate the evolution from BKT to BCS superfluidity 
through the 2D-3D crossover. 

This work is supported by the Australian Research 
Council Centre of Excellence for Quantum- Atom Optics. 
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